Enhancement of superhorizon scale inflationary curvature perturbations 
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We show that there exists a simple mechanism which can enhance the amplitude of curvature 
perturbations on superhorizon scales during inflation, relative to their amplitude at horizon crossing. 
The enhancement may occur even in a single-field inflaton model, and occurs if the quantity a<f)/H 
becomes sufficiently small, as compared to its value at horizon crossing, for some time interval during 
inflation. We give a criterion for this enhancement in general single-field inflation models. 



PACS numbers: 
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I. INTRODUCTION 

The standard, single-field, slow-roll inflation model 
predicts that the curvature perturbation on comoving 
hypersurfaces, 1Z C , remains constant from soon after the 
scale crosses the Hubble horizon, giving the formula 



TZ C « K c (t k ) 



2nd 



(1) 



iH 



where H is the Hubble parameter, d is the time deriva- 
tive of the inflaton field d, and t k is a time shortly after 
horizon crossing. However, one may consider a model 
in which slow-roll is violated during inflation. Recently, 
Leach & Liddle studied the behavior of the curvature 
perturbation in a model in which inflation is temporarily 
suspended, finding a large amplification of the curvature 
perturbation relative to its value at horizon crossing for a 
range of scales extending significantly beyond the Hubble 
horizon. 

In this short paper, we consider single- field inflation 
models and analyze the general behavior of the curva- 
ture perturbation on superhorizon scales. We show an- 
alytically when and how this large enhancement occurs. 
We find that a necessary condition is that the quantity 
z = ad/ H becomes smaller than its value at the time 
of horizon crossing. We then present a couple of inte- 
grals which involve the above quantity and which give a 
criterion for enhancement. 



II. ENHANCEMENT OF THE CURVATURE 
PERTURBATION 

We assume a background metric of the form 



ds' 



= -dt 2 + a 2 (t)S i jdx i dx j 
= a 2 (r)) (-di] 2 + 6,,,l.v'<l.i J ) 



(2) 



On this background the growing mode solution of the cur- 
vature perturbation on comoving hypersurfaces is known 



OU-TAP-155, PU-RCG 01/04, |astro-ph/0101406| 



to stay constant in time on superhorizon scales in the ab- 
sence of any entropy perturbation |l|,00-0] . This follows 
from the equation for 1Z C 



2—Tl' + k 2 K c = 0, 

z 



(3) 



where the prime denotes the conformal time derivative, 
d/dr], and z — ad/H. One readily sees that on superhori- 
zon scales, when the last term can be neglected, there 
exists a solution with 7Z C constant, which corresponds to 
the growing adiabatic mode. 

However, this does not necessarily mean that 1Z C must 
stay constant in time after its scale crosses the Hubble 
horizon. In fact, if the contribution of the other inde- 
pendent mode (i.e. the decaying mode) to 1Z C is large 
at horizon crossing, 1Z C will not become constant until 
the decaying mode dies out. The important point here 
is that the decaying mode is, by definition, the mode 
that decays asymptotically in the future, but it does not 
necessarily start to decay right after horizon crossing. In 
what follows, we show that there indeed exists a situation 
in which the decaying mode can stay almost constant for 
a while after the horizon crossing before it starts to de- 
cay. In such a case, the contribution of the two modes 
to the curvature perturbation is found to almost cancel 
at horizon crossing. This gives a small initial amplitude 
of 1Z C , but results in a large final amplitude for TZ C after 
the decaying mode becomes negligible. 

Let u{rj) be a solution of Eq. (Q) for any given k. For 
much of the following discussion it is not necessary to 
specify the nature of the solution u, but for clarity let 
us identify it straightaway as the late-time asymptotic 
solution at fy* (taking 77* for instance as the end of in- 
flation). For any other solution, v(rj), independent of 
u(rf), it is easy to show from Eq. (g) that the Wronskian 
W = v'u — u'v obeys 



-2-W. 

z 



w 

and hence W oc 1/z 2 . Therefore we have 



(4) 



1 



,vy = W 



(5) 



Hence the decaying mode, v, which vanishes as i] — > rj* 
may be expressed in terms of the growing mode, it, as 



v(rj) oc u(rj) 



drj' 



n. z 2 {rj')u 2 {jj') 



Without loss of generality, we may assume v = u at some 
initial epoch, which we take to be shortly after horizon 
crossing, rj — rjk (<?/*)■ Then v is expressed as 



v(rj) = u(rj) 



D(vk) ' 



where 



DM 



Jri 



z 2 (rj')u 2 (rj') ' 



(6) 



(7) 



and, for convenience, the conformal Hubble parameter 
7~Lk = (&'/a)fc at rj = rj k is inserted to make D dimen- 
sionless. In terms of u and v, the general solution of 1Z C 
may be expressed as 



U c (rj) = au(rj) + 0v(rj), 



(8) 



where a and (3 are constants and we assume a + /3 = 1 
without loss of generality. Thus, if the amplitude of 1Z C 
at horizon crossing differs significantly from that of the 
growing mode, au(rj k ), it can only be because \(3\ 3> 1. 

Using Eq. (||) and noting a + /3 = 1, H c and 1Z' C at the 
initial epoch rj = rju are given by 



K' c (Vk) = u'( Vk ) 



3(1 - a)H k u(r] k ) 
D k 



(9) 



where Dj~ = D(rjk). Then a can be expressed in terms of 
the initial conditions as 



1 + D k 



1 



3H k 



K 



u 
u 



(10) 



If we assume TZ c {qk) to be a complex amplitude deter- 
mined by an initial vacuum state for quantum fluctu- 
ations, then 1Z' C / '{TLkTlc) at the time of horizon cross- 
ing will be at most of order unity. This implies that 
| a |, and hence \j3\, can become large if ^ 1 or 
{D k /H k )\u'/u\ » 1. 



III. LONG- WAVELENGTH APPROXIMATION 

Equation (|j|) can be written in terms of the canonical 
field perturbation, Q — zTZ, as 



k A Q = 0. 

z 



(11) 



From this we see that the general solution for k 2 -C \z"/z\ 
is given approximately by 



TZ C w A + B 



di 



(12) 



where A and B are constants. 

The requirement that v — > as rj — > ry, uniquely identi- 
fies the decaying mode as proportional to f n * dr/' /z 2 (i]') 
in Eq. (|l^), but one is always free to include arbitrary 
contributions from the decaying mode in the growing 
mode. Nonetheless, it is convenient to identify the con- 
stant A in Eq. ( |l2| ) as an approximate solution for the 
growing mode, u, on sufficiently large scales. Thus we 
put the lowest order solutions for u and v as 



uq = const. 



V = U Q 



(13) 



where and in the rest of the paper D (rj) is the integral 
given by Eq. (Q) but with u approximated by uq, 



Jr) 



Z 2 (j)>) 



(14) 



As long as the slow-roll condition is satisfied, the above 
approximate solutions are accurate enough. However, the 
next order correction to the growing mode u may become 
substantial if there is an epoch at which the slow-roll 
condition is violated JtJ. It then becomes important to 
make clear how one chooses the growing mode. 

Rewriting Eq. (|^) in an iterative form, with the lowest 
order solution u — uq, a growing mode solution is given 
by 



u = [1 + F(ri)} u 
F(rj) = k 2 



di rv 



„ z 2 ( V >) 



*V)— (15) 
u 



where the boundary condition is chosen so that u — > u Q 
as rj — > rj*. The above equation tells us when and how 
the approximation u « uq on superhorizon scales may be 
invalidated once the 0(k 2 ) effect is taken into account. 

In the long-wavelength approximation, F can be ap- 
proximated by setting u = uq in the integral to obtain 
the solution to 0(k 2 ) accuracy. Thus we take the ap- 
proximation, 



F( v ) 



dr)' 



z 2 {ri")dri" 



(16) 



The 0(k 2 ) effect cannot be neglected if this integral be- 
comes larger than unity. As may be guessed from the 
form of the integral, such a situation appears if there is 
an epoch during which z 2 {rf) -C z 2 {rj k ). 

To be specific, let us assume z{rj) <C z k — z(rj k ) for 
i] > r)o(> rj k ). Then F{rj) will become large and approx- 
imately constant for rj k < rj < rjo and will decay when 



2 



t] > r) . Incidentally, this behavior is quite similar to the 
behavior of the lowest order decaying mode vo(rj) given 
in Eq. (|T^). In other words, the growing mode can be 
substantially contaminated by a component that behaves 
like the decaying mode, and it can no longer be assumed 
as being constant on large scales. 

The above discussion suggests that we may take ad- 
vantage of the ambiguity in defining the growing mode 
to redefine u accurate to 0(k 2 ) by 



= [1 + F(r,)} u -F k v (v), 



(17) 



where F k = F(jjk)- The growing mode will be now ap- 
proximately constant on super horizon scales: u ~ uq, or 
at least u(j]k) = u(r)*). However, u'/u at 77 = r/k will no 
longer be negligible. We find 



-F k 



Di 



-F k . 



Then Eq. @ for a may be approximated as 



1 



3Tik Tic 



F k , 



(18) 



(19) 



where D k and Fk are those given in the long-wavelength 
approximation, Eqs. ([l4]) and (^), and for definiteness 
we will take (k/H k ) 2 = 0.1. 

In slow-roll inflation, the time variation of is small 
and z increases rapidly, approximately proportional to 
the scale factor a. Hence neither the integral Dk nor 
Fk cannot become large. Soon after horizon crossing 
TZ' C /TZ C <C H, so that a w 1 and the standard result 
TZ c (r]) w 7Z c (r]k) holds. However, if the slow-roll condi- 
tion is violated, may become very small and z may 
decrease substantially to give a large value of Dk and 
F).. (The case where z actually crosses zero is treated 
separately in an Appendix.) Then at late times, we have 



Tl c (v*) = au(?7*) ~ ctu(r/k) = aTZ c (r]k) ■ 



(20) 



Thus the final amplitude will be enhanced by a factor 
\a\, which can be large if Dk 3> 1 or Fk 3> 1. 



IV. STAROBINSKY'S MODEL 

As an example we consider the model discussed by 
Starobinsky J|] , where the potential has a sudden change 
in its slope at = 4>o such that 



_ , V + A + (<j> - <f> ) 

v ^)-\ Vo + A _^_ (j)o) 



for 
for 



> 00 
< 00 



(21) 



If the change in the slope is sufficiently abrupt j£| then 
the slow-roll can be violated and for A + > A_ > the 
field enters a friction-dominated transient (or "fast-roll" ) 
solution with sa —3770 || until the slow-roll conditions 
are once again satisfied 



3H 1 



For < 



-A + for > 0o 

-A- - (A+ - A_)e- 3H ° A * for < 0o 



we have 



-a - 



A e H„At + ( A+ _ A _y-2H At 



3iJ 2 



(22) 



(23) 



This decreases rapidly to a minimum value z m i n ~ 
(A-/A + ) 2 / 3 zq for A + 3> A—, which can cause a signifi- 
cant change in TZ C on superhorizon scales. 

First let us discuss the behavior of D(rf). For a mode 
that leaves the horizon in the slow-roll regime z grows 
proportional to a while > 0o, so that the integrand of 
D{rj) remains small. Hence D{rf) w Dk, which implies 
TZ c (j]) TZ c {rjk) until r\ = 770- Even after the slow-roll 
condition is violated lZ c {v) still remains constant until 
z becomes smaller than Zk and the integrand of -D(?y) 
becomes large again. Then D(r\) may increase rapidly 
until Tl c approaches the asymptotic value for 77 — > 77*, 
given by Eq. (J2fj). Substituting the above solution for z 
in Eq. (||) into Eq. @ we obtain 



Di 



(24) 

which shows that for A+/A- 3> 1, we have Dk > 1 on 
scales (A-fA+y^Ho <k< {A+fA^^Ho. 

A similar behavior is expected for F(rj). Using again 
the solution for z in Eq. ([23]) , the double integral in 
Eq. (|l5|) is evaluated to give 
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for k > (k/H k )H 
for k < {k/H k )Ho 



(25) 



Thus F k > 1 for (A^/A+Y^Ho <k< (A + /A-)^ 5 H . 

Combining the effects of Dk and F k , we see that the 
correction due to Fk dominates on scales k < Ho and Dk 
on scales k > Ho ■ In particular the spiky dip in the spec- 
trum seen in Fig. at k ~ (A-IA+f^Ho is caused by 
Fk, i.e., it is the 0(k 2 ) effect in the perturbation equa- 
tion (|3|). To summarize, the curvature perturbation is 
significantly affected by the discontinuity at ~ 0o even 
on superhorizon scales from k ~ (A- / A + ) x ^ 2 Ho up to 
k ~ (A+fA-Y^Ho- 

Similar behavior was observed in the model studied by 
Leach & Liddle || for false- vacuum inflation with a quar- 
tic self- interaction potential || , whose power spectrum is 
shown in Fig. |^. In this model there is no discontinuity 
in the potential, so the oscillations seen in Starobinsky's 
model are washed out. 
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FIG. 1. The power spectrum for the Starobinsky model M 
with A+/A- = 10 4 . Plotted are the exact asymptotic value 
of the curvature perturbation TZ^.(r/ r ), the horizon-crossing 
value IZciVk), an d the enhanced horizon-crossing amplitude 
o?"Rj^{r)k) using the long-wavelength approximation. The 
range of scales between the dotted lines corresponds to modes 
leaving the horizon during the transient epoch, defined as the 
region where z' j z < 0. Also plotted is the slow- roll amplitude 
K* given by Eq. (pi). 



In both cases our analytic estimate of the enhancement 
on superhorizon scales is in excellent agreement with the 
numerical results on all scales. Thus our approximate 
formula for a given by Eq. ( |l9| ) will be very useful for es- 
timation of the curvature perturbation spectrum in gen- 
eral models of single-field inflation. 

It may be noted that in the Leach & Liddle model the 
long-wavelength condition, fc 2 <C \z" / z\, is violated for 
modes fc < 7^o- It is rather surprising that our long- 
wavelength approximation still works very well for this 
model. 



V. INVARIANT SPECTRA 

A striking feature of these results is that the modes 
which leave the horizon during the transient regime share 
the same underlying spectrum as that produced during 
the subsequent slow-roll era. This is a manifestation of 
the 'duality invariance' of perturbation spectra produced 
in apparently different inflationary scenarios |lC| ]. 

Starting from a particular asymptotic background so- 
lution, z(rj), one finds a two parameter family of solutions 



z{rj) = Cxzirj) + C 2 z(r)) 



Z 2 (l]') ' 



(26) 



which leave z" / z unchanged in the perturbation equa- 
tion (^) and thus generate the same perturbation spec- 
trum from vacuum fluctuations flQj (up to the over- 
all normalization C±). The variable z itself obeys the 
second-order equation 
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FIG. 2. Power-spectrum for the false-vacuum quartic 
model as in Fig. |lj 



z" + a 



,d 2 V 



5H 2 + W + 2 



n 



12 



'U 2 



z = 0. (27) 



Thus for a weakly interacting field (d 2 V/d(j) 2 « constant) 
in a quasi-de Sitter background (H w constant) the equa- 
tion can be approximated by the linear equation of mo- 
tion 



z"+la 



,d 2 V 



2H' 



0. 



(28) 



The general solution z(r]) is related to the asymp- 
totic late-time solution z(rf) by the expression given in 
Eq. @. 

This means that the usual slow-roll result [taking 
w —(dV/d<f>)/3H] for the amplitude of the curvature 
perturbations in Eq. (|l|) 



3ff 3 \ 
27r(dV/#) J 



(29) 



may continue to be a useful approximation even when the 
actual background scalar field solution at horizon cross- 
ing is no longer described by slow-roll, as was noted pre- 
viously by Seto, Yokoyama and Kodama [p"lf and seen in 
our figures. 



VI. SUMMARY 

In summary, we have studied the enhancement of the 
curvature perturbation on superhorizon scales possible in 
some models of inflation. We have found that the curva- 
ture perturbation can be enhanced on superhorizon scales 
even in single-field inflation, provided that the slow-roll 
condition is violated and a<j>/ H becomes small compared 
to its value at horizon crossing. We have presented a 
quantitative criterion for this enhancement, namely that 
either of the integrals Dk and Fk defined by Eqs. (0) 
and (O) , respectively, becomes larger than unity. In the 



long- wavelength approximation (k 2 -c \z" / z\) these in- 
tegrals are expressed in terms of the background quan- 
tity z = a<fi/H, as given by Eqs. (|lj and (|l6|), so an 
analytical formula for the final curvature perturbation 
amplitude may be derived without assuming slow-roll in- 
flation. In the case of a weakly self-interacting field in 
de Sitter inflation we recover the usual slow-roll formula 
for the amplitude of the scalar perturbations even when 
the background solution is far from slow-roll at horizon 
crossing. 
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for 77 < rjQ. Thus introducing the function D{rf) by 



D{rj) 
3Hk 



lim 



dij 



i^k^k 



2, ,2 



(hj 



i z k u k 



Vo+e 

2 zM 



(A7) 



and Dk = D(rjk), where uq — u(rjo), the decaying mode 
v normalized to u at 77 = rjk is given by 



v[r}) = u(r?)-~— 



(A8) 
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Thus exactly the same argument applies to this case, by 
replacing the original Dk by the above Dk ■ 



APPENDIX A: IF </> CROSSES ZERO 

The case when cf> and hence z changes its sign can be 
treated as follows. For simplicity, let us assume z changes 
the sign only once at 77 = 770- Since the integral Fk is still 
well-defined in this case, we focus on the integral Dk- 



In the vicinity of 77 = 770, 
z'o(V ~ Vo) where z' = z'(r/ ) 
becomes 

" d 2 2d 

drf 



z can be expressed as z = 
Hence the equation for 1Z C 



TZ r = 0. 



77 - 770 drj 

The two independent solutions can be found as 



C I 



D 



1 



v - m 



Vo) 



(Al) 

(A2) 
(A3) 



It is apparent that u should be chosen as the growing 
mode, and it remains constant across the epoch 77 = 770- 
We require v to describe the decaying mode. As before, 
we consider an integral expression of v in terms of z 2 and 
u. Then 



drj 



_drf_ 

z 2 C 2 



(A4) 



for 77 > 770. This v behaves in the limit 77 

1 

v 



Vo 



as 



^(77-770) 



(A5) 



This should be extended to the region 77 < 770 as the 
solution (A3), which implies 



u lim 



dr] 1 



+ 



dry' 



rio+e 



7 -C 2 e 



(A6) 
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